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Abstract 

> 

\/~) . For eleven examples of one-dimensional quantum mechanics with shape-invariant 

potentials, the Darboux-Crum transformations in terms of multiple pseudo virtual state 
wavefunctions are shown to be equivalent to Krein-Adler transformations deleting mul- 
tiple eigenstates with shifted parameters. These are based upon infinitely many poly- 
nomial Wronskian identities of classical orthogonal polynomials, i.e. the Hermite, 
Laguerre and Jacobi polynomials, which constitute the main part of the eigenfunctions 
of various quantum mechanical systems with shape-invariant potentials. 
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*u '. 1 Introduction 



The virtual state wavefunctions are essential for the construction of the multi-indexed La- 
guerre and Jacobi polynomials [U [2] . They are polynomial type solutions of one- dimensional 
Schrodinger equations for shape-invariant potentials [31 IU [5]. They are characterised as 
having negative energies (the groundstate has zero energy), no zeros in the physical domain 
and that they and their reciprocals are square non-integrable. By dropping the condition of 
no zeros and the reciprocals are required to be square-integrable, pseudo virtual state wave- 
functions are obtained. In most cases, the virtual and pseudo virtual state wavefunctions 
are obtained from the eigenfunctions by twisting the parameter(s) based on the discrete 
symmetries of the Hamiltonian [Tj. Starting from a shape-invariant potential, a Darboux 



transformation [6j [7] in terms of a pseudo virtual state wavefunction (p{x) with energy 8 
produces a solvable system with an extra eigenstate below the original groundstate with 
energy £ and eigenfunction y5(x) -1 . This method of generating a solvable system by "adding 
an eigenstate" below the groundstate is known for many years, starting from the simplest 
harmonic oscillator potential examples [8] and followed by many authors [9] [H]. As re- 
marked by Adler (15] for the harmonic oscillator case and generalised by the present authors 
[T6] for other potentials, such a system can be derived by special types of Krein-Adler trans- 
formations. That is, the Krein-Adler transformation for a system with negatively shifted 
parameters in which the created state will be the groundstate. The transformation use all 
the eigenstates between the new and the original groundstates. 

In this paper we present straightforward generalisation of the above result for various 
shape-invariant potentials listed in section [3j Coulomb potential with the centrifugal barrier 
(C), Kepler problem in spherical space (K), Morse potential (M), soliton potential (s), Rosen- 
Morse potential (RM), hyperbolic symmetric top II (list), Kepler problem in hyperbolic 
space (Kh), hyperbolic Darboux-Poschl- Teller potential (hDPT), on top of the well-known 
harmonic oscillator (H), the radial oscillator (L) and the Darboux-Poschl- Teller potential 
(J). They are divided into two groups according to the eigenfunction patterns in § 13 . 11 We 
mainly follow Infeld-Hull jl] for the naming of potentials. A Darboux-Crum transformation 
in terms of multiple pseudo virtual state wavefunctions is equivalent to a certain Krein- 
Adler transformation deleting multiple eigenstates with shifted parameters. In contrast to 
the use of genuine virtual state wavefunctions PQ, not all choices of the multiple pseudo 
virtual states would generate singularity free systems. The singularity free conditions of the 
obtained system are supplied by the known ones for the Krein-Adler transformations |15j . 

Underlying the above equivalence are infinitely many polynomial Wronskian identities 
relating Wronskians of polynomials with twisted parameters to those of shifted parameters. 
These identities imply the equality of the deformed potentials with the twisted and shifted 
parameters. This in turn guarantees the equivalence of all the other eigenstate wavefunctions. 
We present the polynomial Wronskian identities for Group A; the harmonic oscillator (H), 
the radial oscillator (L) and the Darboux-Poschl- Teller potential (J) and some others. For 
Group B, the identities take slightly different forms; determinants of various polynomials with 
twisted and shifted parameters. The infinitely many polynomial Wronskian identities are the 
consequences of the fundamental Wronskian (determinant) identity (I2.12p as demonstrated 
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in section HI 

This paper is organised as follows. The essence of Darboux-Crum transformations for the 
Schrodinger equation in one dimension is recapitulated in § 12.11 The definitions of virtual 
states and pseudo virtual states are given in § 12.21 In section [3] two groups of eigenfunction 
patterns are introduced in § 13. II and related Wronskian expressions are explored in § 13.21 The 
details of the eleven examples of shape-invariant systems are provided in § l3.3h § l3~l~3l Section 
H] is the main part of the paper. We demonstrate the equivalence of the Darboux-Crum 
transformations in terms of multiple pseudo virtual states to Krein-Adler transformations in 
terms multiple eigenstates with shifted parameters. The underlying polynomial Wronskian 
identities are proven with their more general determinant identities. The final section is for 
a summary and comments. 

2 Darboux-Crum & Krein-Adler Transformations 

Darboux transformations in general [5] apply to generic second order differential equations 
of Schrodinger form 

d 2 

U = -— + U(x), U^(x)=8^(x) (S,U{x)eC), (2.1) 
ax 1 

without further structures of quantum mechanics, e.g. the boundary conditions, self-adjointness 
of "H, Hilbert space, etc. In the next subsection, we summarise the formulas of multiple Dar- 
boux transformations, which are purely algebraic. 

2.1 General structure 

Let {<fj(x), £j} (j = 1,2, .. . , M) be distinct solutions of the original Schrodinger equation 

H<pj(x) = E m {x) (£j 6 C ; j = 1, 2, . . . , M), (2.2) 

to be called seed solutions. By picking up one of the above seed solutions, say <fi(x), we 
form new functions with the above solution ip(x) and the rest of {(p k (x), £ k } (k ^ 1): 

,[!]/ v def W[(^i,-0](x) _ ^{x)d x j}{x) - dxlf^X^jx) [i] def W[^i,^ fc ](x) 

V\{x) (P!(X) (p!(x) 

It is elementary to show that ip^(x), ip^ 1 (x) (= f <^i(x) _1 ) and (p^\(x) are solutions of a new 
Schrodinger equation of a deformed Hamiltonian %W 

n [1] = + U [1] (x), U [1] (x) = U(x) - 2^1og|^(x)|, (2.4) 



with the same energies £, £\ and 

U [1] i> [1] (x) = £^{x), H [1] ^\x) = £ x ^ l {x), (2.5) 

rt\f k {x) = £ h <pW{x) (Ml)- (2.6) 

By repeating the above Darboux transformation M-times, we obtain new functions 

m def W^^.-.^M.fflx) (2 7) 

W[(p 1 ,(p2,...,<PM](x) 



^ g WK»...,ft,..,»,l(,) ( J = l,2,..., M), (2.* 



which satisfy an M-th deformed Schrodinger equation with the energies £ and £j [7J [15] : 



= -£. + uM(x), U^(x) ^U(x)-2dl\og\W[ Vl ^ 2 ,...,<p M }(x)l (2.9) 

n [M] ^ [M] (x)=£^ M \x), n [M] iff I] (x)=£ I] (x) (j = l,2,...,M). (2.10) 

Here W[/i, f 2 , • • • , /nK^) is a Wronskian 

'd^Ux)- 



W[/ 1 ,/ 2 ,...,/ n ](x) d ^ f det(^^) 



l<j',fc<n 



For n = 0, we set W[-](x) = 1 and W[/i, f 2 , ■ ■ ■ , fj, ■ ■ ■ , f n ](x) means that fj(x) is excluded 
from the Wronskian. In deriving the determinant formulas fl2.7p - fl2.8l) and (12.91) use is made 
of the properties of the Wronskian 

W[gh, gf 2 , gf n ](x) = g(x) n W\f u / 2 , . . . , f n ](x), (2.11) 
W[W[f 1 J 2 ,...J n ,g],W[f 1 ,f 2 ,...J n ,h]](x) 

= W[A, f 2 , . . . , f n ](x) • W[A, f 2 , . . . , f n , g, h](x) (n > 0). (2.12) 

Another useful property of the Wronskian is that it is invariant when the derivative -4 is 
replaced by an arbitrary 'covariant derivative' Di with an arbitrary smooth function qi(x): 

A = ^ - <a(x), W[A, / 2 , . . . , f n ](x) = det( J D J _ 1 • • • D 2 D x f k {x)\^ n , (2.13) 

with -Dj-i • • • D 2 D\ | ._.,= 1. Under the change of variable x — > 77(2:), the Wronskian behaves 

f j {x)=F J (7i(x)) } WfA,/^...,/^.) = (^) ln(M W[^^,..,^(^)). (2.14) 
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Obviously the potential U^ M \x) (12. 9p is independent of the order of the seed solutions. 
The zeros of the seed solution ip\{x) (the Wronskian W[<fi, . . . , (fj^](x)) induce singularities 
of the potential U^(x) in (U^(x) in q> ). 

We apply Darboux transformations for various extensions of exactly solvable potentials 
U(x) G IR in one dimensional quantum mechanics defined in an interval X\ < x < X2- We 
assume that the potential is smooth in the interval and the system has a finite (or an infinite) 
number of discrete eigenstates with a vanishing groundstate energy: 

H(f>„{x) = £ n (p n {x) (0 < n < n max or n e Z> ), (2.15) 
= £ < £i < £2 < • ■ ■ , 

dx (f) m (x)4> n (x) = h n 5 mn , h n >0. (2.16) 



^mi 0n) 



xi 



Hereafter we consider real solutions only and use the term eigenstates in their strict sense, 
i. e. they only apply to those with square integrable wavefunctions and correspondingly the 
eigenvalues. The zero groundstate energy condition can always be achieved by adjusting the 
constant part of the potential U(x). Then the Hamiltonian is positive semi-definite and it 
has a simple factorised form expressed in terms of the groundstate wavefunction 4>o(x) which 
has no node in (xi, 22): 

H = + U(x) = A* A, U(x) = ^Pf^- = (d x w(x)) 2 + d 2 x w{x), (2.17) 

ax z <Po{ x ) 

A def d _ QatW , x s A \ = _d _ d xW f x \ w ( x \ G R Q ( X ) = e -(-). (2.I8) 
dx dx 

For all the examples of solvable potentials to be discussed in this paper, the main part of the 
discrete eigenf unctions {4> n (x)} are polynomials P n (r)(x)), in a certain function r)(x), which 
is called the sinusoidal coordinate [17]. The original systems to be extended usually contain 
some parameter(s), A = (Ai,A2, ...) and the parameter dependence is denoted by "H(A), 
^4(A), £„(A), 4> n (x; A), P n (rj(x); A) etc. We consider the shape-invariant [3] original systems 
only, which are characterised by the condition: 

A(X)A(\) ] = A(X + <5)U(A + S) + £i(A), (2.19) 

or equivalently 

[d x w(x; A)) 2 — dlw(x; A) = [d x w(x; A + S)) 2 + d 2 x w(x;\ + S)+£i(X). (2.20) 
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Here S is the shift of the parameters. 

The explicit extensions depend on the choices of the seed solutions {(fj(x),Sj} (j = 
1, . . . , M). The obvious choices are a subset of the discrete eigenfunctions {(f)j(x), £j} (j G V), 
in which T> is a subset of the index set of the total discrete eigenfunctions. By using those 
from the groundstate on O , lv . . , 4>m-i (T> = {0,1, . . . , M — 1}) successively, Crum [7] has 
derived an essentially iso-spectral extension 

H [M] <j)l M] (x) = £ n ^\x) (n = M,M + l,...), (2.21) 



^ ^■■■■■^y;) , e 1 ) = n & - ^) • (2.22) 

W[0 O ,01, • • • 

C/M( x ) d ^ f U(x) - 2d 2 x log|W[0 o , 0i, ••• , <Pm-i}(x)\, (2.23) 

in which the potential C/W (x) (s = 1, 2, . . . , M) at each step is non-singular. Shape-invariance 
means simply 

U [s \x; A) = U(x; X + sd) + £ S (X) (s = 1,2, . . . , M), (2.24) 
<P [ : ] +n {x; A) oc n (a;; A + sS) (s = 1, 2, . . . , M; n = 0,l,...,). (2.25) 

By allowing gaps in the final spectrum, Krein and Adler (15] have generalised Crum's results 
for V = f {d\, d 2 , . . . , cIm} {dj G Z> : mutually distinct): 

rt M ^\x) = E n <ff\x) {niV), (2.26) 

e^)- W ^' 0d ;-"-'^'t; ]( W (4 M UL M] ) = f[(^-^)-M m „, (2.27) 

tfM^) = C/(x) - 2«9, 2 log|W[0 dl , <p d2 , <P dM ](x)\. (2.28) 



M-l 



The potential U^ M \x) is non-singular when the set T>, which specifies the gaps, 0^ (x) = 
(dj G T>), satisfies the conditions [TBI H~6]: 

M 

Krein-Adler conditions : [J(m - dj) > ( Vm G Z> ). (2.29) 

These simply mean that the gaps are even numbers of consecutive levels. These extensions are 
well-known. In the next subsection, we consider seed functions which are not eigenfunctions. 
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2.2 Virtual Sz pseudo virtual states 

Now let us consider extensions in terms of non-eigen seed functions. For simplicity sake, 
we first consider the cases of exactly iso-spectral extensions (deformations). The seed func- 
tions {<fj(x),Sj} (j = 1,2,..., M) satisfying the following conditions are called virtual state 
wavefunctions: 

1. No zeros in x\ < x < X2, i.e. <fj(x) > or <fj(x) < in X\ < x < x 2 . 

2. Negative energy, Sj < 0. 

3. <fj(x) is also a polynomial type solution. 

4. Square non-integrability, (<fj,<fj) = oo. 

5. Reciprocal square non-integrability, (cpj ,cpj ) = oo. 

Of course these conditions are not totally independent. The negative energy condition is 
necessary for the positivity of the norm as seen from the norm formula (12.271) . since a similar 
formula is valid for the virtual state wavefunction cases when Sd 3 is replaced by Sj. 

When the first condition is dropped and the reciprocal is required to be square integrable 
(ipj , (pj ) < oo, such seed functions are called pseudo virtual state wavefunctions. When 
the system is extended in terms of a pseudo virtual state wavefunction (fij(x), the new 
Hamiltonian TV- 1 ' has an extra eigenstate Lp^ l (x) with the eigenvalue £j, if the new potential 
is non-singular. The extra state is below the original groundstate and "H^ is no longer 
iso-spectral with %. This is a consequence of (I2.5p . Its non-singularity is not guaranteed, 
either. When extended in terms of M pseudo virtual state wavefunctions {(pj(x),£j} (j = 
1,2, .. . ,M), the resulting Hamiltonian V\ M ^ has M additional eigenstates fij(x) ( 12. 8ft . if 
the potential (x) is non-singular. They are all below the original groundstate. 

Since tpj(x) is finite in X\ < x < x 2 , the non-square integrability can only be caused by 
the boundaries. Thus the virtual state wavefunctions belong to either of the following type 
I and II and the pseudo virtual state wavefunctions belong to type III : 



Type I : 




(2.30) 



7 



rxi+e i-X2 

Type II : / dx(pj(x) 2 = oo, / dx <Pj(x) 2 < oo, 

rxi+e PX2 

& / dxtpj(x)~ 2 < oo, / dxcpj(x)~ 2 = oo, (2.31] 

Ai J X2~ e 



rxi+e fX2 

Type III : / dx(fj(x) 2 — oo or / <ix <^-(x) 2 = oo, 

>/ xi J X2~- e 

pxi+e PX2 

& / dx (pj(x^) 2 < oo, / dx (fj(x)~ 2 < oo. (2.32) 

</iei J X2~- e 

An appropriate modification is needed when X2 = +00 and/or xi = —00. Hereafter we denote 
the pseudo state wavefunctions by {0 V }? which are the main ingredients of this paper. 

The Darboux-Crum transformations in terms of type I and II virtual state wavefunctions 
have been applied to achieve exactly iso-spectral deformations of the radial oscillator poten- 
tial and the Darboux-Poschl- Teller potentials [I], which generate the multi-indexed Laguerre 
and Jacobi polynomials. 

In this paper we show that the Darboux-Crum transformations in terms of a set of 
V = f {di, d 2 , . . . , d M } pseudo virtual state wavefunctions 

WW^ft). (2-33) 
= W wS'^'"' , ^t ]( f ) ' (eUf^n^n-^)-^, (2.34) 

W[0 dl ,0 d2 ,...,0rf M J(x) j=1 

d 4 f [/( x ) - 2 «9 2 log|Wfo, d2 , . . . , dM ](x)|, (2.35) 

are equivalent to a system generated by Krein-Adler transformations f l2.26p - fl2.28p specified 
by a certain set T> ( 14. 2 p of eigenfunctions with shifted parameters for various systems with 
shape-invariant potentials, which are listed in the subsequent section. In contrast to the 
extensions in terms of the virtual states, these extensions in terms of pseudo virtual states 
are not iso-spectral and the obtained systems are not shape-invariant. 

As seen in each example, the virtual and pseudo virtual wavefunctions are generated 
from the eigenstate wavefunctions through twisting of parameters based on the discrete 
symmetries of the original Hamiltonian. 
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3 Examples of Shape- invariant Quantum Mechanical 
Systems 



Here we provide the essence of shape-invariant and thus exactly solvable one- dimensional 
quantum mechanical systems. The first five examples § l3.3f § l3~7n have infinitely many discrete 
eigenstates, whereas the rest § l3.8f § l3~T3l has finitely many eigenstates. They are divided into 
two groups of eigenfunction patterns. The eigenfunctions {4> n } and the pseudo virtual state 
wavefunctions {0 V } 

H(X)<f) n (x; A) = £ n (\)(j) n (x] A) (0 < n < n max (A) or n 6 Z> ), (3.1) 
H{ A)0 v (x; A) = £ v (\)Mx; A) (v G V( A)) , (3.2) 

are listed in detail for reference purposes. Here V(A) is the index set of the pseudo virtual 
state wavefunctions, which is specified for each example. The type I and II virtual state 
wavefunctions for systems with finitely many eigenstates will be discussed in a separate 
paper [IB] . 

The basic tools of shape-invariant systems (I2.19p -f l2.20p are the forward shift and back- 
ward shift relations: 

A{X)(j> n {x] A) = / n (A)0 n _ 1 (ar; A + 8), (3.3) 
^(A)Vn-i(x; A + 8) = 6 n _ 1 (A)0 n (z; A), (3.4) 
£ n (X) = f n (\)b n ^{\). (3.5) 

For the parameters of shape-invariant transformation A — > A + 8, we use the symbol g 
to denote an increasing (8 = 1) parameter and the symbol h for a decreasing (8 = — 1) 
parameter and \x for an unchanging [8 = 0) parameter, except for the Darboux-Poschl- 
Teller potentials (J) in which g and h are both increasing parameters. Throughout this 
paper we assume that the parameters g and/or h take generic values, that is not integers or 
half odd integers. 

In the next subsection, we introduce two groups of the eigenfunction patterns of the eleven 
examples of shape-invariant systems, § l3.3h § l3~T3"l Then the basics of the Wronskians for the 
two groups are discussed in § 13.21 The explicit formulas of the eigenvalues, eigenfunctions, 
pseudo virtual state wavefunctions and the discrete symmetries of the eleven shape-invariant 
systems are provided in § l3.3h § l3TT3l for reference purposes. Among them, we also report type 
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II virtual state wavefunctions for two potentials (C) § 13.61 and (Kh) § 13.121 These have been 
reported in [llj and [H], respectively in connection with the M = 1 rational extensions. We 
stress that multi-indexed orthogonal polynomials can be constructed for these two potentials 
in exactly the same way as in |lj. For another type of virtual state wavefunctions for the 
potentials with finitely many discrete eigenstates § l3.8f § l3~T3| see a subsequent publication 
[18J. The sections |3~3TI3. 131 could be skipped in the first reading. 



3.1 Two groups of eigenf unctions patterns 



In all the eleven Hamiltonian systems § l3.3h § l3~T3l the energy formula of the pseudo virtual 
state is related to that of the 'eigenstate' with a negative 'degree': 

£ V (X)=£^ 1 (X). (3.6) 

These Hamiltonian systems are divided into two groups of eigenfunction patterns: 

po(x; A) : Group A 



^>Jx: A) 



X)P n ( v (x);X) : 



Vn)0&; A) 



(3.7) 



po(x; A + nS) : Group B 

Seven potentials (H), (L), (J), (M), (s), (hst) and (hDPT) belong to group A and four 
potentials (C), (K), (RM) and (Kh) belong to group B. The group A has a simple structure. 
The ! n' dependence is carried only by the degree of the polynomials (except for (M)). They all 
satisfy the closure relations p2] and rj(x) is called a sinusoidal coordinate. The groundstate 
wavefunctions of this group satisfy 



<> (x; X + S) 
(j) (x; A) 



dx 



1 


H, s, hst 


2 


L 


-4 


J 


-1 


M 


4 


hDPT 



(3.f 



The group B has a more complicated structure. The 'n' dependence appears also in the 
other factor and in the parameter of the polynomials a n , (3 n . 

The pseudo virtual state wavefunctions are obtained from the eigenstate wavefunctions 
by twisting A — >■ t(A) (and x — > ix for (H) and (L)): 



K(x; A) 



<Po(v)(x;X) 



v (x;t(A)) = 4>o( v )(x; X)^ v (r](x); A), 

Po(x; t(A)) : Group A _ , , 

P (x;t(X)+vd) : Group B ' A) - t(A)j : 



(3.9) 
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with a slight modification for (H) and (L): 

(H) : v (x; A) = i~ v (f) v (ix) = 4>o( v )(x)^ v (r](x)) , 

0o(v)(V) = <fio(ix), £v0?) = i~ v P v (ir]), (3.10) 
(L) : v (a;;A) = i*~V v (iac;t(A)) = 0o(v)(^; A)£ v (r/(a;); A), 

0o(v)(z;A) =^-Vo(^;t(A)), ^(^A) = P v (-r / ;t(A)). (3.11) 

The twist operation t satisfies t(A + ad) = t(A) — aS (a G C). 

Corresponding to the forward and backward shift relations of the eigenfunctions ( 13. 31) - 
(13.4j) and the energy factorisation formula ( 13.51) . those for the pseudo virtual state wavefunc- 
tion read: 

A{\)Mx;\) = -eb v {-\)i v+1 {x;\ + 6), e = | ~* " , (3.12) 

A(X)^ Y+1 (x; A + 6) = - e 'f v+1 (-\)Mx; A), e' = j j , (3.13) 

£ V (A) = ee / / v+1 (-A)6 v (-A) = /_ v _x(A)6_ v _ 2 (A) = £_ V _ X (A). (3.14) 

3.2 Wronskian formulas 

Based on the eigenfunction patterns, the Wronskians of the eigenfunctions and the pseudo 
virtual state wavefunctions for the set T> = {di, d 2 , . . . , g?m} 

W[(j> dl ,(f) d2 , (j>d M ]{x; A), W[0 dl , (p d2 , <t> dM ](x; A) 

are reduced to determinant formulas of the polynomials. Since the latter is obtained from 
the former by twisting, we present derivations of the former. 

For Group A, the reduction to the Wronskian of the polynomials is achieved by the 
Wronskian formulas (l2TTTj) and ^Wf (P n (x; A) = P n (r](x); A)) : 

W[0 dl , fa, . . . , (j> dM ]{x- A) = A) M W[P dl , P d27 . . . , P dM }(x; A) 
/dn(r) \ hM(M-l) 

= o (x; A) M (^) W[P dl ,P d2 , . . . ,P dM ](v(x); A). (3.15) 

For Group B, the derivative operator d^~ x in the Wronskian is replaced by 'covariant 
derivatives' A{\ + (j - 2)6) ■ ■ ■ A(X + S)A(X) (12TT31 and use is made of the forward shift 
relation (13.31) to obtain 



J-2 

.A(A + (j - 2)<5) ■■■A(\+ S)A(X)4> n (x- A) = J] / n _i(A + ^6) ■ </> n - j+ i(x; A + (j - 1)<5) 

i=0 
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J-2 

= II /«-*( A + ^) ' ^ A + (?7^); A + (j - 1)5). (3.16) 

i=0 

We have 



l<j,k<M 



W[(f) dl ,(f) d2 ,...,(f) dM }(x] A) = det(jJ/ dfc _;(A + i5) • (f) dk - j+1 (x; A + (j - 1)5)) 

j=0 

M j-2 

IT 0o(^; A + 4^) • det (TT f dk ^{\ + ^5) • P dk _ j+1 (n{x); A + (j - 1)5)) . (3.17) 



fc=l i=0 



This method is applicable to Group A, too, in which YlkLi 4>o{ x '-, A + 45) is replaced by 

Jj0o(a:;A + O-l)A) = o (z; A) M (^/-^J 
j'=i 

Let us summarise the results: 

W[0 dl ,<f> d2 ,..., 4> dM ] (x; A) = Axi(x; A)Sx>(?7(x); A), (3.18) 



^(x;A) d ^ f <| ^ - > — (3 ig) 

rife=i 0o(^; A + 45) : Group B 

- fa A) d = f / C ' tf(tf " llw [ p *- P *' ■ ■ ■ ' P ^](^ A) : Group A 
Idet^A))^.^ : GroupB ' 

A) = J] + ^) • P dk - j+ i{v, A + (j - 1)5). (3.21) 

j=0 

The pseudo virtual state wavefunction Wronskian is simply obtained by the twisting: 

W[4> dl , 4> d2 , 4> dM }(x; A) = A v (x; A)S 27 ( ? 7( a ;); A), (3.22) 

A v (x;X)^\ h^^) M (^r MiM - 1] :Grou P A ^ 
' Ilfcli 0o (a:; t(A) + 45) : Group B 



[x . x) M {c -i^mM-D ;GroupA 



~ X) d 4 f I C ^ M(M 1)w ^*' • • • ' ^/](^; A) : Group A 24 
IdetlX^A))^^ : Group B ' 



j-2 



A ) = II /*-* W*) + i*) ■ U-i+ifa A - 0' - 1)5). (3.25) 



i=0 

(Note that the expressions of E v and for Group B are valid for Group A, too.) Their 
relations are 

Av(x; A) = A v (x;t(X)), E V ( V ;\) = E v ( V ;t(X)), (3.26) 

12 



with a slight modification for (H) and (L): 

H : A v (x; A) = A v (ix), E v (t]; A) = i~ t7> E v (iri) , 

L: A v (x;X)=i^^ M A v (ix;t(X)), Ed(t]; A) = z M(M-1) Hu(-ry; t(A)). 
Both Hp (77; A) and Hp (77; A) are polynomials in r] and their degrees are generically £t>'- 

M 

^ = J>--M(M-1). (3 - 27) 

Remark : Strictly speaking, the notation T> of the polynomials Hp and Hp represents an 
ordered set. By changing the order of d/s, Hp and Hp may change sign. On the other hand 
the functions A-p and Ad are invariant under the permutations of d/s. In order to avoid 
excessive appearance of ± signs in the general formulas involving Hp and Hp, etc, we adopt 
the following convention. The formulas in §|4] involving the Wronskians and determinants of 
various polynomials depending on V, V and other sets are understood to be true up to a 
multiplicative constant ±1 coming from the multi-linearity of the determinants. This does 
not affect the main Propositions, e.g. (I4.19p - fl4.20p in Proposition 14.21 

3.3 Harmonic oscillator (H) 

The well known system of the harmonic oscillator has infinitely many eigenstates: 

A : none, 8 : none, —00 < x < 00, 

w(x\ A) = — -x 2 , U(x; A) = x 2 — 1, 

£ n (A) = 2n, 7](x) = x, f n (X) = 2n, & n -i(A) = 1, 

4> n (x;X) = (f) (x;X)P n (ri(x);X), <j) (x; A) = e"^ 2 , P n (r]; A) = H n {rj) , 

h n (X) = 2 n n\^. 

Here H n (x) is the Hermite polynomial. The system satisfies the closure relation introduced in 
[T7] . The pseudo virtual state wavefunctions are obtained by the following discrete symmetry: 

n(x) = -H(x)\ x _ Hx +s_ 1 (x), 

V (>; A) = i' v (f) v (ix; A) = e^ x2 i~ v H v (ix) (v G Z> ), 
S V (X) = -£ V (A) + £^{X) = £_ v _i(A). 
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3.4 Radial oscillator (L) 

The radial oscillator potential has also infinitely many discrete eigenstates in the specified 
parameter range: 

\ = g, 6 = 1, < x < oo, g > -, 

w(x; A) = -\x 2 + g\ogx, U(x; A) = x 2 + 9 ^ ~ l) - (1 + 2g), 
2 x z 

S n {X) = 4n, v {x)=x 2 , f n {X) = -2, 6 n _ 1 (A) = -2n, 

4>n{x\ A) = 4>o(x; X)P n (r](x)] A), o (x; A) = e"^ 2 x 9 , P n (r/; A) = L ( n~^{v), 

h n (X) = ^Y{n + g + \). 

Here L^a\ri) is the Laguerre polynomial. The system satisfies the closure relation [T7| . 
There are three types of discrete symmetries: 

type I : U{\) = -H(\)\ x ^ ix -2(1 + 2#), 

type II: H{\) = H(t(\)) + 2(1 - 2g), t{\) = l-g, 

type III : U{\) = -n(t(X))\ x ^.+£^{X), 

and the pseudo virtual states are generated by type III: 

v (x; A) = z'-^iz; t(A)) = e^V^-r^); t(A)) (v e Z> ), 
4(A) = -£ v (t(A)) +£_x(A) = £_ V _!(A). 

The type I and II virtual states are obtained by using type I and II discrete symmetries pQ. 

3.5 Darboux-P6schl- Teller (DPT) potential (J) 

The DPT potential has also infinitely many discrete eigenstates in the specified parameter 
range: 

71 3 

\={g,h), S =(1,1), 0<x<-, g,h>-, 

/ x\ ,, tt, x\ gig-l) h(h-l) , , . o 

A) = g log smi + fi log cos 2, (7 (x; A) = 5 1 5 (g + h) , 

sin £ cos" 1 x 

£ n (A) = 4n(n + g + /i), rj(x) = cos 2x, / n (A) = — 2(n + # + /i), &n-i(A) = — 2n, 
n (x; A) = o (x; A)P n (r?(x); A), o (x; A) = (sinx^cosxf, P n (ri; A) = P^~ iM) (7?), 
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hJX) 



T{n + g + \)T{n + h+\) 



2 n\ (2n + g + /i)r(n + g + h) ' 
Here P^'^irj) is the Jacobi polynomial. The system satisfies the closure relation [T7j . 
There are three types of discrete symmetries: 

type I: H{\) = H(i\\)) + (1 + 2g){l - 2h), t\\) = (g, 1 - h), 
type II : H{\) = H{l l \X)) + (1 - 2g){l + 2h), t n (A) = (1 - «/, h), 
type III: "H(A) = "H(t(A)) + £_i(A), t=po£, t(A) = (1 - g, 1 - h), 

and the pseudo virtual states are generated by type III: 

v (x; A) = v (:r; t(A)) = (sin x^cosx) 1 -^ (r}(x); t(A)) (0 < v < g + /i - 1), 
£ V (A) = £ v (t(A)) +£_ X (A) = £_ v _i(A). 

The type I and II virtual states are obtained by using type I and II discrete symmetries pQ. 
The Hamiltonian % has also the 'left-right' mirror symmetry x — >• f — x, g h. 

3.6 Coulomb potential plus the centrifugal barrier (C) 

The system has also infinitely many discrete eigenstates in the specified parameter range: 

A = g, 6 = 1, 0<x<oo, 9>^, 

»(*;A)= 9 log*-H U(x ,X) = ^l- 2 - + \, 
g x A x g z 

E n {\) = \ - - — 7](x) = x" 1 , / n (A) = — — -^—r, 6„_i(A x - 



g 2 (g + n) 2 '' //(// • " 2g 

n (x;A) =e-l^x° +n P n {r ] (x) ] \), o (x; A) = e"f x\ P n ( V] A) = V n L^ {^V' 1 ), 

k nW =(^y 9+ ^n2g + n) . 

The discrete symmetry and the pseudo virtual state wavefunctions are: 
H(\) = H(i(\))+£- 1 (\), t(X) = l-g, 

Mx; A) = v (x; t(A)) = e^x^^P^x); t(A)) (0 < v < g - 1), 
£ v {\) = £ v (t(A)) + f_!(A) = f-v-i(A). 

For g — 1 < v < 2g — 1 (<? > |), the discrete symmetry generates the type II virtual states 
v (x;A) HI]. 

This system can be obtained from the Kepler problem in spherical space § 13. 71 in a certain 
limiting procedure. 
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3.7 Kepler problem in spherical space (K) 

The system has also infinitely many discrete eigenstates in the specified parameter range: 

3 

\=(g,n), 6 = (1,0), 0<a;<7r, g > -, /i > 0, 



w(x; A) = g log sin x — —x, U (x; A) = ^ 2 ^ — 2/i cot x + ^ — g 2 , 



2 

(7 sin x g 

2 



^n(A) = {g + nf -g 2 + ^- - ^ r](x)=cotx, 

g 2 (g + n) 1 

f (\\ _ ff 2 (ff + ^) 2 + ^ 2 , ^ _ w(2y + n) 
g{g + n) 2 g 

(f) n (x; A) = e~^ x (siax) 9 n P n (rj(x); A), (f) (x; A) = e~^ x (smx) 9 , 

P n ( m \)=C n Pl a ^Xiri), a n = -g-n+-^—i, n = -g - n - -^—i, 

g + n g + n 

e -^2 1 - 2( -e +n h(g + n)T(2g + n) 



K{X) 



^ ((9 + nf + j^)T{g + ^i)V(g - ^i) 



The discrete symmetry and the pseudo virtual state wavef unctions are: 
H(X) = U{K\))+£-i(\), t(A) = (l-g,ti), 

v (x; A) = v (x; t(A)) = e^ x (smx) l ~ 9+v P v (r](x); t(A)) (0 < v < 2g - 1), 
£ V (X) = £ v (t(A)) + £_i(A) = f_ v _i(A). 

The Coulomb potential plus the centrifugal barrier § 13.61 is obtained by the following 
limit: 

x = lx c ) » = ± hm(f) 2 ^(A)=^ c (A c ), 
together with the eigenfunctions. 



3.8 Morse potential (M) 

The system has finitely many discrete eigenstates < n < n mSuX (X) = [h]' in the specified 
parameter range ([a]' denotes the greatest integer not exceeding and not equal to a): 

A=(/i,/i), S = (—1,0), — oo < x < oo, h,fx>0, 

w(x; A) = hx - fie x , U(x; A) = ^i 2 e 2x - (i(2h + l)e x + h 2 , 

n — 2h 

S n {X) = h 2 -{h- n) 2 , V (x) = e~ x , f n (X) = 6 n _!(A) = -2n/x, 

n (x; A) = o (x; X)P n ( V (x); A), o (x; A) = e hx ~^\ 
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P„(„; A) = (2 w -r»4— >(W>, MA) - egg^j - 

The system satisfies the closure relation [T7] . 

The discrete symmetry and the pseudo virtual state wavefunctions are: 

H{\) = n(t(X)) + £_i(A), t(A) = (-1 - h, -//), 

v (x; A) = v (x; t(A)) = e^+^+^Pv^); t(A)) (v G Z> ), 

f v (A) = £ v (t(A)) + £_!(A) = f_ v _!(A). 

This system can be obtained from the hyperbolic Darboux-Poschl- Teller potential § 13.131 
in a certain limiting procedure. 

3.9 Soliton potential (s) 

The system has finitely many discrete eigenstates < n < n majX (\) = [h]' in the specified 
parameter range: 

A = h, S = —1, — oo < x < oo, h > 0, 

h(h + 1) 

w(x; A) = —h log cosh x, U(x; A) = 5 \~ ^ 2 > 

cosh x 

2 /. ^2 x /xx l t /xn n(2h-n) 



8 n (X) = h - (h-ny, rj(x) = sinhx, /„(A) = h, & n -i(A) 
n (a;;A) = o (:r; \)P n (i](x); A), </> (z;A) = (coshx 



/1 



.-ft 



P^x); A) = (coshx)"p(^ h -")(tanhx), ^(A) = -— 1 + J 

n\ (n — n)L {2h — n + 1) 

The system satisfies the closure relation [T7]. One can rewrite P n {v, A) as 

(/l — [■ !: M-])rii+ii (■_/,_! 

p n 77; a) = ; ' ;; [ 2 1 z"pi »• 

(n — n+ gjp+i] 

where [a] denotes the greatest integer not exceeding a. 

The discrete symmetry and the pseudo virtual state wavefunctions are: 

n(\) = n(i(x))+£- 1 (\), t(A) = -i-/i, 

v (x; A) = v (x; t(A)) = (cosh x) h+1 P v (77(0;); t(A)) (v G Z> ), 
£ V (A) = £ v (t(A)) + f_i(A) = £_ v _i(A). 

This system can be obtained from the Rosen- Morse potential § 13.101 by taking fi — > 
limit. 
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3.10 Rosen-Morse potential (RM) 

This potential is also called Rosen-Morse II potential. The system has finitely many discrete 
eigenstates < n < n max (A) — [h — ^/JI}' in the specified parameter range: 

A=(/i,/i), d = (—1,0), — oo < x < oo, h > y/JI > 0, 

/i h(h + 1) 2 /i 2 

wfx; A) = — /i log cosh x — — x, U(x;X) = 5 h 2/i tanh x + h +717, 

n. cosh a; ft 2 



/!~ /I 2 

£ n (A) = /i 2 - (h - n) 2 + — - — — , r)(x) = tanhrc, 

n z (h — ny 



h 2 (h-n) 2 -fi 2 n(2h-n 

Jn{A) — TTT ^ j On-l{<*) 



h(h — n) 2 ' h 

<fi n (x;\) = e~~^ x (coshx)~ h+n P n (ri(x); A), <fi (x;X) = e~^ x {cosh.x)~ h , 

P n (r ]] X) = P^\r ] ), an = h-n + T ^-, (3 n = h - n 

K{\) 



h — n h — n' 

2 2h - 2 "(h - n)T(h + ^ + l)T(h + 



n\((h-nf-jj^)r(2h-n+l) 

By taking the limit fi — > 0, the soliton potential § 13.91 is obtained. Based on the symmetry 
x — > —x, ji — » — ji, positive \i is selected. 

The discrete symmetry and the pseudo virtual state wavef unctions are: 

U{\) = U(l{\)) + 5_i (A), t(A) = (-1 - h,fi), 

<fi v (x; A) = <f> v (x;i{\)) = e^ a; (coshx) /l+1+v P v (r/(x); t(A)) (v e Z> ), 
£ V (X) = £ v (i(\)) + S^(X) = £_ V _!(A). 

3.11 Hyperbolic symmetric top II (hst) 

The system has finitely many discrete eigenstates < n < n max (A) = [h]' in the specified 
parameter range: 

\=(h,/i), S = (—1,0), — 00 < x < 00, h,fi>0, 

w(x; A) = — h log cosh x — /1 tan -1 sinh x, 

TT( -/i(/i + l)+/i 2 + / u(2/i+l)sinhx 2 

(7 (x; A) — — rj h a , 

cosh a; 

n — 2h 

£ n {\) = h 2 - (h- n) 2 , T](x) = sinhx, / n (A) = — - — , &„_i(A) = -2n, 
<P n (x- A) = o (x; A)P ft (7j(aO; A), A ) = e^ tan " sinhl '(cosha;)- h , 
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P n (ri;X) = r n P^\ir]), a = -h - \ - in, = -h-± + i(i, 

h (x) = 7rT(2h-n + l) 

nl ' 2 2h n\(h-n)T(h-n + \ + iij)T(h-n + \-iny 

The system satisfies the closure relation [T7] . 

The discrete symmetry and the pseudo virtual state wavefunctions are: 

U{\) = U(i{\)) + £_x(A), t(A) = (-1 - h, -n), 

Mx;X) = (j> v (x;t(\)) = e^ tan " lsinh:E (cosha;) /i+1 P v (7 ? (a;);t(A)) (v G Z> ), 
S V (X) = £ v (t(A)) + £_ X (A) = 5_v_i(A). 

3.12 Kepler problem in hyperbolic space (Kh) 

This potential is also called Eckart potential. It has finitely many discrete eigenstates < 
n < n max (A) = [y/jl — g)' in the specified parameter range: 

A=0,/i), S = (1,0), 0<x<oo, s/ji>g>-, 

w(x; A) = g log sinh * - H x , U(x; A) = - 2„ coth x + g 2 + 

g sinh x g 

u 2 a 2 

S n (X) = g 2 - (g + n) 2 + — - - — ■ — — , r)(x) = cothx, 

g 2 (g + n) 1 

f r\\ - V 2 -9 2 (g + n) 2 , /.x _ n(2g + n) 
g(g + n) 2 g 

4> n (x;X) = e~~9+™ x (smh.x) 9+n P n (r](x); A) , 0o(^; A) = e~^ x (smhx) g , 

P n ( m \) = Pi a ^( V ), a n = -g-n+-^— 1 fi n = - g - n - 

g + n g + n 

h (g + n)T(l-g + ^)T(2g + n) 

^ +2n ^(^-(9 + nr)T(g + ^- n y 

The discrete symmetry and the pseudo virtual state wavefunctions are: 

U{\) = U(l{\)) + f_i(A), t(A) = (1 - g,fi), 

v (x;A) = <f) v (x;t{\)) =e^ x (sinhx)^ +v+1 P v (r/(x);t(A)), 

4(A) = £ v (t(A)) + £_ X (A) = f-v-iCA) (0 < v < (/ - 1, v > ^ + <? - 1). 

For g — 1 < v < 2g — 1 > |), the discrete symmetry generates the type II virtual states 
v (x;A) p. 
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3.13 Hyperbolic Darboux-P6schl- Teller potential (hDPT) 

It has finitely many discrete eigenstates < n < n max (A) = [-^Y in the specified parameter 



range: 



*={g,h), S = (1,-1), < £ < oo, h>g>-, 



w(x; A) = glogsinhx — h log cosh x, U(x; A) = ^-^ 2 — + (h — g) 2 

sinh x cosh x 

S n (X) =4n(h- g-n), t)(x) = cosh2a;, f n (\) = 2(n + g-h), & n -i(A) = -2n, 
n (x;A) = o (x; A)P n (r](x); A), O (^;A) = (sinh x) 9 (cosh x)~ h , 

P„(„; A) = Pi'"--*' (,), M A) - r( " + 9 + " 9 " " + ^ 



2n\ (h- g -2n)T{h-n + \)' 

The eigenvalues can be also expressed as £ n (A) = 4((^=^) 2 — — n) 2 )- The system 
satisfies the closure relation [T?] . 

Three types of discrete symmetries are: 

type I : H{\) = H{i\X)) - (1 + 2g){l + 2/i), t r (A) = (g, -1 - h), 
type II : H(A) = H(t n (A)) - (1 - 2<j)(l - 2h), t n (A) = (1 - g, h), 
type III: H(X) = U(i(\)) + £-i(A), t = t n o t 1 , t(A) = (1 - -1 - /i). 

The pseudo virtual state wavefunctions are generated by type III: 

4> v (x; A) = <f) v (x; t(A)) = (sinhx^-^coslixf+^^x); t(A)) (v G Z> ), 
£ V (A) = £ v (t(A)) + f_ x (A) = f_ v _ x (A). 

The type I and II virtual states are obtained by using type I and II discrete symmetries. 
The Morse potential § 13.81 is obtained by the following limit: 

x = ±(x M -a), g=\ne a + a, h = \^e a + 2h M + a, lim H{\) = 4^ M (A M ), 

together with the eigenfunctions. 

4 Main Results 



dcf 



Let us introduce appropriate symbols and notation for stating the main results. Let D 
{d±, d,2, . . . , djvf} (dj £ ^>o) be a set of distinct non-negative integers. We introduce an 
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integer N and fix it to be not less than the maximum of T>: 

N > max(P). (4.1) 

Let us define a set of distinct non-negative integers V = {0, 1, . . . , iV}\{(ii, d 2 , ■ ■ ■ , (1m} to- 
gether with the shifted parameters A: 

V = {0, 1, . . . , di, . . . , d 2 , . . . , d~M, ...,N} = {ei, e 2 , . . . , eAr+i_ A/ }, 

dj = N-dj, X = A- (N + 1)5. (4.2) 

Starting from the well-defined original system ( 13. ip . the system after Darboux-Crum trans- 
formations in terms of a set of pseudo virtual state wavefunctions T> is described by the 
Hamiltonian H DC 



A 2 

U DC (x) = U(x; A) - 2d 2 x logl Wfc , fa, . . . , fa, ](x; A) I . (4.3) 



General theory presented in §[2] states that, if the Hamiltonian H DC is non-singular, the 
eigenstates are given by and $° c : 



dc /_\ W [<M > fa > ■ ■ ■ > fa/ , 0n] (a:; A) 



I T Til (n = 0,l,...,n max (A)oroo), 
W [0 dl ,<p d2 ,. . ., <p dM \ {x; A) 

^pc (a;) = Wfa,^,- : -,^,- : -,^](x;A) (J = 1)2; __ M); 

W[(j) dl ,(/) d2 , . . . ,faj(x; A) 
H DC $° c (x) = £ n (A)$° c (:r), U GC $f c (x) = ^.^(A)^ 00 ^). (4.4) 

The system after Krein-Adler transformations in terms of f> with shifted parameters A is 
described by the Hamiltonian "H KA 

n KA = -^- 2 +u KA (x), 

dx 2 

U KA (x) = U(x; A) - 2dl log|W[0 o , 0i, • • • , fa, • • • , fa,, • • • , Mfa A) | . (4.5) 

Here we assume that the original system (13. ip with the shifted parameters A has square 
integrable eigenstates, etc. If the Krein-Adler conditions are fulfilled, eigenstates are given 
by and $^ A : 

^ {x) = W[0o, 1; ■ ■ ■ , fa, • • • , fa,, ■ --An, ^ +1+ n](x; A) (n = 0j Xj _ jfw(A) OT ^ 
W[0 , 0i, . . . , fa, . . . , fa,, . . . , 4>n]{x; A) 
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£ KA , x = W[0p, 1; . . . , fc, . . . , . . . , d - M , . . . , <j) N }(x; A) . = i 2 

W[0 , 0i, • • • , fen ■ ■ ■ > , • • • , 0jv] (s; A) 
^ KA $^ A (x) = £ N+1+n (\)^ A (x), n KA $f A (x) = £j.(\)<tf A (x). (4.6) 

Note that n max (A) in § 13.84 — § f3TT3l satisfies 

n max (A) - (N + 1) = n max (A). (4.7) 

The differential equations (14. 4p and (14.61) (with n G Z>o) hold irrespective of the non- 
singularity of T-L DC and T-L KA . 

Our main results hold for any one of the shape-invariant quantum mechanical systems in 

Proposition 4.1 The two systems with 'H DC and 1-L KA are equivalent. To be more specific, 
the equality of the potentials and the eigenf unctions read: 

U DC (x)-£_ N _ 1 (\) = U KA (x), (4.8) 
$° c (x) oc $« A (x) (n = 0, 1, . . . ,n max (A) or oo), (4.9) 
$° c (x) cc <lf A (x) (j = l,2,...,M). (4.10) 

The singularity free conditions of the potential are 

N+l-M 

JJ (m-ej)>0 (VmeZ> ). (4.11) 

3=1 

For M — 1, T> — {di}, T> = {0, 1, . . . , di, . . . , N}, the above conditions are satisfied by 
even d±, d\ G 2Z>o In other words, the pseudo virtual state wavefunctions {<p v } for even 
v are nodeless. The above equalities ( up to multiplicative factors ) f 14 . 8 [) - (I4.10p ( (14. 9 p with 
n G l*>o) are algebraic and they hold irrespective of the non-singularity conditions (14.1 ip . 

The following energy formulas satisfied by the eleven systems 

£n(A) — £_7v_i(A) = £ N+ i +n (X), 
£-?-i(\) - f-jv-iCA) = ^at-v(A), (4.12) 

and the illustration (Figure 1) would be helpful to understand the Proposition. The form 
of the energy curve is that of the DPT potential but the situation is similar in all other ten 
potentials. 
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N - dj N + 1 



n 



Figure 1: The left represents the Darboux-Crum transformations in terms of pseudo virtual 
states. The right corresponds to the Krein-Adler transformations in terms of eigenstates 
with shifted parameters. 



If the equality of the potentials (14. 8 p is shown, then, under the condition (14. lip , the 
Hamiltonians "H DC — £_jv-i(A) = "H KA are non-singular and it implies (14. 9 p and ( I4.10p . 
Therefore the remaining task is to show (14. 8p . 

The basic ingredients of the Proposition 14.11 are the two Wronskians 



which determine the potentials (14. 8p and eventually all the eigenfunctions (I4.9p - (l4.10p . By 
using the Wronskian formulas (I3.22p and ( 13. 18j) they are expressed by the determinants of 
polynomials: 



These two polynomials Ex>(rj; X) and S^(r/;A) have the same degree £x> = if) (I3.27p . An 
important relation between Ad(x; A) and Af,(x; A) is that its ratio is independent of M and 
dj in each of the eleven potentials § l3.3h § l3TT3~| which is shown by straightforward calculation. 
We set the ratio as F(x, N, A): 





W[0 dl ,0 d2 , . . .,(f) dM )(x; A) = A v (x; \)E v (r](x); A), 

W[0 O , 0i, • • • , fe> • • • , 0J M , • • • , 0v]O; A) = Af>(x; X)Ef,(r](x); A). 



(4.13) 
(4.14) 



A v (x; A) 
Af,{x- A) 



F(x,N,X). 



(4.15) 
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From this independence and the obvious relation ^4{}(x; A) = A{y(x; A) = 1, the ratio F can 
be expressed as 

F(x, N, A) = .1;,,., v;!.r: A) = -= 1 . -r . (4.16) 

^{o,i,...,jv}(^; A) 

It is easy to show the following identity for any of the eleven systems: 

U(x; A) — 2d 2 x \og\F(x, N, A)| - £_iv_i(A) = C/(x;A). (4.17) 

By this identity, the equality of the potentials (14. 8[) reduces to 

d>g|S P (77(x);A)| ^loglS-p^); A)|. (4.18) 

Since both polynomials are of the same degree £%>, this means that the two polynomials are 
proportional. We state this as 

Proposition 4.2 Two polynomials characterising "H DC and "H KA , Hx?(r/; A) and S^(?7;A) ; 
are proportional: 

SpfaAJocSsfaA). (4.19) 
In particular, it means polynomial Wronskian identities 

W[U , a , • • • , fa A) oc W[P , Pi, • • • , P Sl , ■ ■ ■ , h M ,---,Pn] fa A) , (4.20) 

for all the systems in Group A. Recall that £ v fa A) = P v fa t(A)) (13. 9p (with a slight mod- 
ification for (H) and (L), fl3TT0l) - fl3~TT]) l For the simplest case of M = 1, JV = d x = t (=> 
V = {£}, V = {1, 2, . . . , £}), the Wronskian identities are reported as (A. 22) in [16] for 
the three types of the classical orthogonal polynomials, the Hermite (H), Laguerre (L) and 
Jacobi (J). To the best of our knowledge, the general polynomial Wronskian identities (14.201) 
have not been reported before. 

The proof of Proposition 14.21 is done by induction in M. 
first step : In the first step we prove ( 14. 19|) for M — 1, N > d\ = v: 

£ v faA) (77; A). (4.21) 

Recall that the differential equations (I4.4p and (14. 6p hold, and 

4> v {x; A) = A {Y} (x; A^fax); A), (4.22) 
W[0 O , 0i, , 4>n](x; A) = A { o,i,...,jv}(^; A)H { o,i,...,jv}(^; A), 
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W[0 O , 0i, . . . , 0v, • • • , <Pn](x; A) = A {0 ,i,...^...,jv}O; ^{o,i,...,^...,N}(v, (4-23) 

By substituting (1Q2]) into the differential equation % DC $f c (x) = £_ v _i(A)$? c (x) it 
becomes 

^/ + 2 , ; d x f+ 1 - U(x; A) + £_ v _i(A) ) / = 0, 4.24 



where / = £ v (^(x); A). By substituting ( 14. 23ft into the differential equation % KK $\ 
£v(A)$f A (x) (14. 6 p and using S{o,i,...,jv}(^(^); ^) — constant (see (14.381) ). we obtain 



x 



d 2 J +2d x \og 



G, 



G 



dgG v 9 X G r .d x G v d x G 



dj + {-^ + - - U(x; A) + S,(X) ) f = 0, (4.25) 



where / = S {0)1 v ..^ ^ N} (rj(x); X), G v = A {0h j^ N} (x; A) and G = A {0jl> ... :N} (x; A). By 
using ( I4.15P and (I4.16p . we find A{ v }(x; X)/G v = F(x,N,X) = 1/G. The above equation 
( I4.25P is rewritten as 

d^pXl , %A {y} (x-,\) x) + S,(X)-2dl log F(x, N, A)) / = 0. (4.26) 

A {v} (x;X) V A {v }(x;A) / 

From ( 14121) and ( 14~T71) . we have 

E/(x; A) - £_ v _i(A) = U(x) A) - <^(A) + 2d 2 JogF(x, N, X). 

Thus / = £ v (j](x); A) and / = H/ 01 ^ N \(j](x); A) satisfy the same differential equation. 
Since both of them are polynomials of degree v in 77, they should be proportional for generic 
parameters [19J. 



second step : Assume that ( 14.191) holds till M (M > 1), we will show that it also holds for 
M + 1. 

Before presenting a general proof, we illustrate the outline by taking the simple case of 
Group A ( I4.20p . We have shown M = 1 case: 

e v (^A)ocW[Po,...,Pv,---,Piv](^A). 

By using the algebraic Wronskian identity (I2.12p . we have 

Wf&.&JfaA) cxW[W[P ,...,P Sl ,...,P N ],W[Po,...,Pd 2 ,...,PN]](vA) 

= ±W[P ,...,P Sl ,...,P S2 ,...,P N ](v;X)-W[P ,P 1 ,...,P N ]( V ;X) 
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<xW[Po,...,P d - 1 ,...,P d - 2 ,...,P N ](v,\), 

where we have used W[P , Pi, . . . , P/v](?7; A) = constant (see (I4.38P ). This is the M = 2 
result. For M = 3, we use the M = 2 results to obtain 

f dl (77; A) • W&fc, (77; A) 

= w[w[a,e*],w[a,U]] M) 

ocW[W[Po,...,4>---^ 2 >--- p Ar],W[P ,...,4 ) --->4 ) --->^v]](r7;A) 

= ±w[p ,...,4,---,4,---,4---^](^ A )- w [ p o,...,4>--->^](^ A ) 

cc W[P , ■ ■ ■ , P 3l , ■ ■ ■ , 4, • • • , Pd 3 , • • • , Piv](ry; A) • £^(77; A). 
This establishes the identities for M — 3. Higher M identities follow in a similar way. 
Let us present a general proof. Eq. fl4.19p is equivalent to 

Wfcxjd2> • • • jdj^j A) W[0p, </)!,..., (f>a M , . . . , (/>jy](x; A) 
A v (x;\) X Av(x;\) 

which implies 

W[fa, d2 , • • • , &iJ(z; ^) oc P(x, iV, A)W[0 o , 0i, . . . , . . . , . . . , 0iv](^; A). 
Assume that (14.19p holds till M (M > 1). By using the Wronskian identity (I2.12p . we obtain 
W[fa,. . . ,0 dM _J(x; A) • W[0 dl , . . . Ad M -xAd M Ad M+x ]{x] A) 

= w [w[fa , . . . , , , wfe , . . . , 4>d M - x , ( x ; A ) 

oc W[F(x, N, A)W[0 o , . . . , fa, . . . , h M -n • • • > • • • > <M> 

P(x, iV, A)W[0 o , • • • , fa, • • • , h M -v • • • ' fc^+i' • • • ' ^1] ( x ; ^) 

= p(x, at, A) 2 w[w[0 o , ---,fa,---, h M ^,- ■ ■ , k M , M, 

W[0 O , ...,fa,..., ^a M -i ' • • • ' ^"m+i > • • • ' <I>n]] (s; A) 

= ±p(x, iv, a) 2 w[0 o , • • • , fa, ■ ■ • , fau-i>' ■ ■ ' <?W • • • ' faii+n • • • > <M(a;; ^) 

x W[0 O , • • • , fa, ■ ■ ■ , 0<J M -i> ■■, ^n](x; A) 
oc F(x, N, A)W[0 o , • • • , fa, ■ ■ ■ , h M+1 , <Pn]{x; A) ■ W[fa,. . . , fa^^x; A). 

This leads to 

W[fa,...,4>d M+1 ](x; A) 
A{d lt ...,d M+1 }(x; A) 
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F ( x , ><)^{o,...,d 1 ,...,d M+1 ,...,N}(x-, A) W[0p, . . . , d -, . . . , 0j M+1 , . . . , 0jy](a;; A) 
A {dl; ... >dM+l} (x; A) ^{o,...,Ji,...,d M +i,-,w}( a; ; 
_ W[0p, . . . , (fa , . . . , 0j M+1 , ■ ■ ■ , 0jy] (g; A) 
^{o,...,Ji,...,d M+ i,.-.,Ar}( :r ; A) 

which means 

5 {di > ...,dM+i}(^; A ) a s {o,...,d 1 ,..., ( i M+1 ,... i 7V}( ? 7; A). 

This concludes the induction proof of (14.191) and the proof of Propositions 14.11 and 14.21 is 
completed. 

In order to obtain the rational forms (the ratios of polynomials) of the eigenfunctions 
1?]) — (I4.10p in Proposition 14. 1\ we need the Wronskian expressions of the numerators: 



W[0 dl , 4> d2 , 4>d M , <p n ] (x; A) = A Vtn (x; \)P Vyn (r}(x); A), 

AvAx . x) = f 0o(o)(x; A) (c-M) W-*>~" : Group A ^ 

I rifcii 0o (z; t(A) + <4<5) • 0o(^; A + n<5) : Group B 
Here Pv,n{v'i A) is a polynomial in r\ and its degree is generically £d + M + n. This leads to 
W[4> dl ,...,4> dM ,<t>n](x\ A) Pt>,n(v(. x )'i A) / O (^;A-M<5) : Group A 



W[4> dl ,...,(j> dM ](x;X) E v (rj(x); A) t o (^;A + n<5) : Group B 
Eq. (13181) leads to 

W[<i) dl ,...,<f) dM ,(f) n \{x;X) P-D,n{r](x); A) / o (a;; A + M6) : Group A 



x 



W[0 dl , . . . ,0d M ](x; A) Sx,(?7(x);A) \ o (:r; A + nS) : Group B 
Pv,n(v') A) == S{ dli ... Afi „}(?7; A) : degree = ^{d lr ..,d M ,n} = tv-M + n. (4.29) 

In the present case we have 

W[0 ei , • • • ,0 ew+1 _ M ,0iv+i+n](z; A) _ ^&,jv+i+n(^(^); A) / o (z; A - M<5) : Group A 



X 



W[0 ei ,...,0 ejv+1 _ M ](x; A) Ef,(r](x); A) \ </> (:r; A + n<5) : Group B 

where the degree of Pt>,N+i+n(v, A) is £5 - (N + 1 - M) + (N + 1 + n) = £ v + M + n. 
Therefore (14.91) implies 

Pv,n{V, A) OC Pf, 

,N+l+n 

(77; A) (nGZ> ). (4.30) 

Similarly (I4.10p implies 

s {*,...A,...,i*r}(^ A ) K ^ ^ A) (j = 1,2,..., M). (4.31) 
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Let us introduce T>± based onD= {d\, d%, . . . , c?Af}: 



V± = {d 1 ±l,d 2 ±l,...,d M ±l}. 
The newly introduced polynomials at n = are related to the old ones: 

Pv,o{m a a + <5), 

Pv,o(v, A) oc 5©+ (77; A + <5). 

Based on these formulas, the following relations follows: 

W[(j) dl ,...,(f) dM ,(j)o\{x;X 



(4.32) 



(4.33) 
(4.34) 



fv(x; A) = f log 



mm dj > 2), 



W[0 dl ,...,0 dM ](x; A) 
(^/^(ar; A)) 2 - d 2 J v (x; A) = (^/ c _(x; A + S)) 2 + d 2 J v _(x; A + 6) + S 1 (X), 
W[0 dl ,...,0 dM) o ](a;; A) 



/■d(x; A) = f log 



W[0 dl ,...,0 dM ](x; A) 
(d x f v {x; A)) 2 - c^(x; A) = (d x f v+ (x; A + <5)) 2 + d 2 J v+ (x; X + 6)+ £ 1 (X), 



(4.35) 



(4.36) 



which have the same forms as the shape-invariance relation f)2.20p but they do not mean 
the shape-invariance. The M = 1 case for (M), (RM) and (Kh) was presented as 'enlarged' 
shape-invariance in [T3l fT4] . 

In the rest of this section we provide the proofs of fl4.33l) - fl4.34l) . For f)4.33p . we note that 
the forward shift relation (13.31) can be written as 



d (j) n (x;\) _ <p n -i(x] X + S) 



dx 4>o(x; A) 



(x;A) 



(4.37) 



By using this, we obtain 



W[(f) dl , . . . ,(f> dM ,(f) ](x] A) 



i (s; A) A1+i W 
^(x; X) 



5di 



M+1/_ 1 XM W 



A) 








• > 


00 


- d 


0rfl 


-dx 


00 



(x; A) 
d 



'>d M 



i (x;A) M+1 (-l) M W 

M 



dx 0o J 
^>di-i{x; X + 6) 



x;X) 



fd M W' 



^d M -i(x;X + S) 



o (x;A) ' O ( X; A) 

!) M II A ( A ) • A ) w [^i-i> ■ ■ ■ > 0^-iK*; A + 5). 
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This can be rewritten as 

-{di,...,rf M ,0}( r /( a O; ^) / i \M TT f ( \\ 

^ {dl -i,...4 M -i } {r]{x); A + d) £± 

and the fact that its right hand side is a constant 



b (x] A)A {dl _ lv .. A/ _i } (x; A + 8) 



A 



{dx,...,d M 



,o}(^; A) 



^o(x; A)A {dl _ lv „ Af _ 1} (x; A + 8) _ 
A{ dl ,...,d M ,o}(x;X) 



can be verified for any of the systems listed in § l3.3f § l37L3l This concludes the proof of 
fl4.33p. We remark that f!4.33p is E{ dl ,...,d M ,o}{r)] A) oc H {dl _i ) ... idAf _i}(77; A + 6), which implies 

^{o,i,...,N}(Vi A) oc H{o,i,...„iyr-i}(f7; A + (5) oc • ■ • oc H {0} (?7; A + iV<5) = constant. (4.38) 
For (14.341) . we note that the forward shift relation for 4> v (x; A) (13.121) can be rewritten as 



d 4> v (x; A) 
dx (po( x 'i A) 



-e& v (-A> 



<) v+1 (x; X + S) 
(j) (x;X) 



(4.39) 



By using this, we obtain 



W[<f> dl ,...,4> du ,<f)o\{x;X) 
= o (x;A) M+1 W 



'0 ' 00 



{x; A) M+1 (-1) M W 



(x; A) 

dx (fio' ' fix 0o 



(x; A) 



M 



x- X) M +\-l) M w\-eb dl {-X) ^+ l( f ; A ;+ 5) , . . . , -e6 dM (-A)^ +l(x; X + 6) 



(x; A) 



io(ar; A) 



(x 



Y[ e fc dj (-A) ■ <f) (x; A)W[0di+i, • • • , 0d M +i] ( x ; ^ + 



M 



This can be rewritten as 

P {du ..., dM} ,o{y{x); X) = yr b 

and the fact that its right hand side is a constant 



o (x; A)A {dl+ i v .. A/+ i } (x; A + (?) 
A{d lt ...,d M },o(x; A) 



f>o(s; A)v4 {dl+lv „ idM+1} (x; A + S) = 

-4{di,...,dM},o(^; A) 



can be verified for any of the systems listed in § l3.3h § l3~l~3l This concludes the proof of 

(EES}. 
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5 Summary and Comments 



In the context of rational extensions of solvable potentials, the concept of the pseudo virtual 
state wavefunctions is introduced. They are obtained by relaxing two conditions of the virtual 
state wavefunctions; the reciprocals are square-integrable and they need not be nodeless. A 
Darboux-Crum transformation in terms of a pseudo virtual state wavefunction will produce 
a new eigenstate below the original ground state. The same system can be derived by a 
special type of Krein-Adler transformation with negatively shifted parameters. The main 
results of the paper is the equivalence of Darboux-Crum transformation in terms of multiple 
pseudo virtual states and Krein-Adler transformations in terms of multiple eigenstates with 
negatively shifted parameters. This is based on polynomial Wronskian identities, which are 
generalisations of those reported by the present authors [TB] a few years ago. The equivalence 
holds for most of the known shape-invariant potentials consisting of eleven explicit examples 
having finite as well as infinite discrete eigenstates. 

The type II virtual state wavefunctions have been obtained by the discrete symmetries 
for two potentials (C) § 13.61 and (Kh) § 13.121 They have been used in the context of M — 1 
rational extensions of these potentials in [TTJ and [H]. Multi-indexed extensions of these two 
potentials can be constructed in exactly the same way as in [T] . In a separate publication [TB] , 
we will discuss rational extensions in terms of genuine virtual state wavefunctions for shape- 
invariant potentials having finitely many discrete eigenstates. They have different features 
from those having infinitely many eigenstates, which have been explored in connection with 
the multi-indexed orthogonal polynomials pQ. 

The multi-indexed Laguerre and Jacobi orthogonal polynomials are labeled by the multi- 
index P, but different multi-index sets may give the same multi-indexed polynomials, e.g. 
eqs.(50)-(51) in [IJ. The proposition 14.21 gives its generalisation. By applying the twist based 
on the type II discrete symmetry to (14.191) . the l.h.s becomes the denominator polynomial 
with multiple type I virtual state deletion and the r.h.s. becomes that of type II. 

After completing the manuscript, we became aware of a recent work [20J, which dis- 
cusses some rational extensions of the harmonic oscillator. They correspond to some special 
examples of the equivalence for the harmonic oscillator (H) and M = 2. 
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